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We include a generalized infinite class of quark-gluon vertex dressing diagrams in a study of how 
dynamics beyond the ladder-rainbow truncation influences the Bethe-Salpeter description of light 
quark pseudoscalar and vector mesons. The diagrammatic specification of the vertex is mapped 
into a corresponding specification of the Bethe-Salpeter kernel, which preserves chiral symmetry. 
This study adopts the algebraic format afforded by the simple interaction kernel used in previous 
work on this topic. The new feature of the present work is that in every diagram summed for the 
vertex and the corresponding Bethe-Salpeter kernel, each quark-gluon vertex is required to be the 
self-consistent vertex solution. We also adopt from previous work the effective accounting for the 
role of the explicitly non-Abelian three gluon coupling in a global manner through one parameter 
determined from recent lattice-QCD data for the vertex. With the more consistent vertex used 
here, the error in ladder-rainbow truncation for vector mesons is never more than 10% as the current 
quark mass is varied from the u/d region to the b region. 

PACS numbers: 12.38.Aw; 11.30.Rd; 12.38.Lg; 12.40.Yx 



I. INTRODUCTION 

In recent years there has been significant progress in 
the study of the spectrum of hadrons, and their non- 
perturbative structure and form factors, through ap- 
proaches that are manifestly covariant and which ac- 
commodate both dynamical chiral symmetry breaking 
(DCSB) and quark confinement Covariance provides 
efficient and unambiguous access to form factors 
Consistency with chiral symmetry and its spontaneous 
breaking is obviously crucial to prevent the pseudoscalars 
from artificially influencing the difficult task of describing 
and modeling the infrared dynamics; this is a role better 
left to other hadronic states that are not so dominated by 
chiral symmetry. The associated concept of a constituent 
quark mass is important and it is often implemented in 
models as a constant mass appearing in the propagator; 
however this idealization runs into trouble for higher ly- 
ing states where the sum of the constituent masses is 
below the hadron mass. This difficulty is marginally evi- 
dent with the p, but it is inescapable by the time one has 
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reached the ground state axial vector mesons (e.g., a%, b\ 
mesons) 

In reality, solutions of the QCD equation of motion for 
the quark propagator (quark Dyson-Schwinger equation 
(DSE)) give a momentum-dependent quark mass func- 
tion. Model calculations, mostly in Landau gauge, typi- 
cally yield a mass function that evolves from the current 
mass value at ultra-violet spacelike momenta to a value 
some 0.4 GeV larger in the deep infrared 0. The propa- 
gator is a gauge-dependent object and the gauge depen- 
dence of this phenomenon has not been fully explored. 
In the chiral limit, such an enhancement is DCSB. At fi- 
nite current mass, models also strongly suggest that the 
enhancement is the same mechanism as DCSB which has 
an important influence over the low-lying hadron spec- 
trum. In the chiral limit, the scalar term of the quark 
self-energy, which shows most of the momentum depen- 
dence, plays a dual role as the dominant invariant am- 
plitude of the chiral pion Bethe-Salpeter equation (BSE) 
amplitude at low momenta Q • In any process where the 
spatial extent of the pion plays an important role, the 
running of the quark mass function is likewise crucial to 
an efficient symmetry-preserving description. Otherwise 
a theoretical model is fighting symmetries. An example is 
provided by the pion charge form factor above the chiral 
symmetry breaking scale, i.e., Q 2 > m p . It is this large 
value of the dressed quark mass function at low spacelike 
momentum that leads, in model solutions of the quark 
DSE, to \p 2 \ M 2 (p 2 ) within a significant domain of 
timelike momenta where these models can be trusted. 
For example, this is sufficient to prevent spurious qq pro- 
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duction thresholds in light quark hadrons below about 
2 GeV 0. 

The task of maintaining manifest covariance, DCSB, 
a running quark mass function and an explicit substruc- 
ture in terms of confined quarks is often met by models 
defined as truncations of the DSEs of QCD 0,11 0. For 
practical reasons the equations must be truncated to de- 
couple arbitrarily high order n-point functions from the 
set of low order n-point functions used to construct ob- 
servables. A common truncation scheme is the rainbow- 
ladder truncation. Here the one-loop gluon dressing of 
the quark (with bare gluon-quark vertices) is used self- 
consistently to generate the quark propagator. In general 
the kernel K of the Bethe-Salpeter equation is given in 
terms of the quark self-energy S by a functional rela- 
tion dictated by chiral symmetry |Toj . This preserves the 
Ward-Takahashi identity for the color singlet axial vector 
vertex and ensures that chiral pseudoscalars will remain 
massless, independent of model details. With a rainbow 
self-energy, this relation yields the ladder BSE kernel. To 
go beyond this level, one needs to realize that the exact 
quark self-energy is given by the same structure except 
that one of the gluon-quark vertices is fully dressed. It is 
the vertex dressing that generates the terms in K beyond 
ladder level. This is the topic we are concerned with in 
this paper. 

The ladder BSE for meson bound states is an integral 
equation with a one-loop kernel structure that must al- 
low for the spinor structure of propagators and the meson 
amplitudes. With the four-dimensional space-time that 
one must use to maintain manifest covariance, and with 
dynamically generated quark propagators that one must 
use to preserve the Ward-Takahashi identities of chiral 
symmetry, the numerical task is large. Any scheme for 
corrections to the ladder truncation will in general add 
the complexity of multiple loop Fcynman diagrams in- 
volving amplitudes that are only known after solution. 
For practical reasons the studies that have been able 
to investigate hadron states be yon d ladder-rainbow (LR) 
truncation in recent years 0, 0, 0, 0] have exploited 
the simplifications following from use of the Munczek- 
Nemirovsky (MN) model [15|- In this case the basic el- 
ement is a delta function that restricts the exchanged 
(or gluon) momentum to zero; it reduces both the quark 
DSE and the meson BSE to algebraic equations. There 
is only one parameter — a strength set by m p . 

This simplified kernel has no support in the ultra-violet 
and one must be wary of its use for related physics. 
Bound state masses are relatively safe in this regard, 
even heavy quark states that sample short distance or 
large momenta due to the large quark mass scale present. 
Even with the MN model, the DSE solutions for the 
quark propagators have the correct power law behavior, 
and support connected to the current quark mass, in the 
ultra-violet, apart from log corrections. The dominant 
qualitative features of DSE solutions of realistic model 
are preserved in the MN model: large infrared strength 
giving DCSB and the (confining) absence of a mass pole. 



Our analysis is not aimed to provide a serious represen- 
tation of experimental data; rather we aim at some un- 
derstanding, even if it is quite qualitative, of the relative 
importance of classes of higher-order diagrams for the 
BSE kernel for bound states. Because of the inherent 
complexity brought by use of a momentum distribution 
as a kernel, there is little information available in the lit- 
erature on this topic. To obtain such information, we feel 
the price paid by dispensing with a clear connection to 
perturbative QCD is worthwhile in the initial stages. 

The first study of the correction to ladder-rainbow 
truncation was made in this context in Ref. where 
a one-gluon exchange dressing of the quark-gluon ver- 
tex was implemented for pseudoscalar and vector mesons 
and scalar and axial vector diquark correlations. Subse- 
quently [H it was realized that the algebraic structure 
allowed a recursive implementation of the ladder series 
of diagrams for the quark-gluon vertex as well as an im- 
plementation of the corresponding series of diagrams for 
the chiral symmetry-preserving BSE kernel. So far as we 
are aware this was the first solution of a BSE equation 
in which the kernel contained the effects from an infinite 
number of loops. In these works the chiral pseudoscalars 
remained massless independent of the model parameter, 
m p received corrections of order 10% from ladder dressing 
of the vertex, and the diquark states evident at ladder- 
rainbow level were removed from the spectrum by the 
dressing effects. The influence of vertex dressing upon 
the quark propagator was also studied. 

There is very little in the way of guidance from realis- 
tic non-perturbative non-Abelian models of the infrared 
structure of the quark-gluon vertex. It has often been 
assumed, e.g., see Ref. that a reasonable beginning 
is the Ball-Chiu 01 or Curtis-Pennington [l^ Abelian 
Ansatz times the appropriate color matrix. These 
Abelian descriptions of the momentum dependence sat- 
isfy the Abelian vector Ward-Takahashi identity, and 
their use makes the implicit assumption that this might 
be a good enough approximation to the corresponding 
identity for QCD, namely the Slavnov- Taylor identity for 
the color octet vertex . The use of an explicit ladder 
sum for the gluon vertex provides easy access to the chi- 
ral symmetry preserving BSE kernel and receives some 
motivation from the fact that a ladder-summed photon- 
fermion vertex combines with the rainbow approxima- 
tion for the fermion propagator to preserve the Ward- 
Takahashi identity for that vertex. 

However when initial results from lattice-QCD simula- 
tions of the gluon-quark vertex became available [20ll2l| . 
it was realized |22| that the color algebra generated 
by any ladder sum for this vertex gives a magnitude 
and strength for the dominant amplitude at zero gluon 
momentum that is qualitatively and quantitatively in- 
compatible with the lattice data and incompatible with 
the leading ultra-violet behavior of the one-loop QCD 
Slavnov- Taylor identity. The infrared vertex model de- 
veloped in Ref. 22] made an extension of the fact that the 
one-loop QCD color structure introduced by the three- 
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gluon coupling repairs the deficiency of a purely ladder 
structure. The color structure of the ladder class of dia- 
grams produces a weak repulsive vertex, while the color 
structure of the three-gluon coupling contribution pro- 
duces an attractive contribution that is enhanced by a 
factor of -iVf at the purely one-loop level. 

These observations from Ref. |22j were blended with 
the algebraic features afforded by the MN model to re- 
examine the relation between vertex dressing, the chi- 
ral symmetry-preserving BSE kernel, and the resulting 
meson spectrum and diquark correlations [T^j . This 
approach introduced one extra parameter (besides the 
gluon 2-point function strength and the quark current 
mass) — an effective net color factor fitted to lattice-QCD 
data on the gluon-quark vertex. The net attraction in 
the vertex, driven by the explicitly non-Abelian 3-gluon 
coupling, had a marked effect: the ladder-rainbow trun- 
cation made m p 30% too high compared to the solution 
from the completely summed vertex. In other words, 
the attraction produced by summed vertex dressing in a 
non-Abelian context is more important than previously 
thought. However in that approach, the structure of the 
vertex is such that the coupling of any internal gluon 
line to a quark, is itself bare. This is not self-consistent 
and one can question what effect this omitted infinite 
sub-class of vertex dressing and BSE kernel contributions 
may have upon the hadron spectrum. 

In the present work we extend the analysis of Ref. 0] 
by the incorporation of a wider class of vertex dressing 
diagrams. We allow the coupling of any internal gluon 
line to a quark to be described by the dressed vertex at an 
order consistent with a given total order in the final ver- 
tex. In the limit of the vertex summed to all orders, this 
becomes the use of the self-consistent quark-gluon vertex 
at every internal location in a diagram. We borrow from 
previous work the use of the MN model of the 2-point 
gluon function to generate an algebraic structure and we 
again incorporate the important non-Abelian three-gluon 
coupling through the device of an effective net color fac- 
tor refitted to the lattice data for the vertex. We use 
the infinite series of diagrams for the BSE kernel gener- 
ated from the chiral symmetry-preserving relation to the 
quark self-energy. We investigate the resulting spectrum 
of pseudoscalar and vector mesons. 

In Section II we describe general properties of the 
quark-gluon vertex and the relationship with the asso- 
ciated BSE kernel that preserves chiral symmetry. Infor- 
mation from the Slavnov- Taylor identity for the gluon- 
quark vertex and the Ward-Takahashi identity for the 
color singlet axial-vector vertex is summarized for rel- 
evance to present considerations. We discuss diagram- 
matic summations that have been used previously to 
model the gluon vertex and the generalized class of di- 
agram considered here. In Section III we introduce the 
interaction model that allows an algebraic analysis and 
we present consequent results for the gluon-quark vertex 
and the self-consistent dressed quark propagator. The 
associated symmetry-preserving BSE kernel is presented 



also. Section IV contains a presentation and analysis of 
the methods and results for the meson masses. In Section 
V there is a summary. 



II. THE QUARK-GLUON VERTEX AND THE 
BETHE-SALPETER KERNEL 

We employ Landau gauge and a Euclidean metric, 
with: {7^,7^} = 2<5 M „; 7+ = 7^; and a ■ b = £) i=1 a % b % . 
The dressed quark-gluon vertex for gluon momentum k 
and quark momentum p can be written ig t c T a (p + k, p) , 
where t c = A c /2 and A c is an SU(3) color matrix. In 
general, T (T (p + k,p) has 12 independent invariant ampli- 
tudes. We are particularly concerned in this work with 
the vertex at k = 0, in which case the general form is 

T a (p) = ai(p 2 )7< T + a 2 (p 2 )7 -pp a - a s {p 2 )ip a 

+ ai(p 2 )i-i a 7 -p (1) 

where on{p 2 ) are invariant amplitudes. In the model 
studies of Refs. [l2j and [LlJ that we build upon, one 
finds CH4 = 0; this will also be the case here. 

As we shall discuss later, we wish to utilize the func- 
tional relation that enables the BSE kernel to be gener- 
ated from the quark self-energy so that chiral symmetry 
is preserved. This requires the vertex to be represented 
in terms of a set of explicit Feynman diagrams. Some 
exact results are known for the vertex at 1-loop order in 
QCD [23. In Landau gauge and to 0(g 2 ), i.e., to 1-loop, 
the amplitude T a is given by 

v£Xp+k,p) = z lla + T^{ P +k,p) + r^ A (p + k, P ), (2) 

with 

T^(p + k,p) = -(C F - ^f) J g 2 D^(p - q)^ 

x S (q + k)^ a So(q)7„, (3) 



and 



T^ A (p + k,p) = -^ J g 2 ^S (p-qh„D^,(q + k) 
xiT% l<r {q + k,q)D v , v {q), (4) 

where — J A d A q/{2i:) i denotes a loop integral reg- 
ularized in a translationally-invariant manner at mass- 
scale A. Here Zi(fj, 2 ,A 2 ) is the vertex renormalization 
constant to ensure T a = r y a at renormalization scale \x. 
The following quantities are bare: the three-gluon ver- 
tex ig / abc Tfjl a (q + k,q), the quark propagator So(j>), 
and the gluon propagator D^(q) = T llu (q)Do(q 2 ), where 
Tuu(q) is the transverse projector. The next order terms 
in Eq. are C(<? 3 ) : the contribution involving the four- 
gluon vertex, and 0{g A ): contributions from crossed-box 
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and two-rung gluon ladder diagrams, and 1-loop dress- 
ing of the triple-gluon vertex, etc. The color factors in 
Eqs. ©and J3J are given by 



t a t b t a = (C F - — )t b = —t b 

2N r 



t ajabc t b = CA u c = El U c 

J 2 2 



t a t a = C F 1 C 



2N n 



(5) 



In contrast, for the color singlet vector vertex, i.e., for 
the strong dressing of the quark-photon vertex, one has 
the one-loop Abelian result 



fW(p + k,p) = Z lla -C F / g 2 D„ v (p-q)^S (q+k) 

Jq 

x la S (qhu- (6) 

To motivate the approximate vertex used in the present 
study, we note that the local color SU(3) gauge invariance 
of the QCD action gives the Slavnov- Taylor identity 0] 
for the gluon vertex 

k^iT^p + k,p) = G(k 2 ) { [1 - B(p, k)] S{ P + k)- 1 

-SW-^l -B(p,k)]}, (7) 

which relates the divergence of the vertex to the quark 
propagator S(p), the dressing function G(k 2 ) of the 
ghost propagator ~G(k 2 )/k 2 , and the ghost-quark scat- 
tering kernel B(j>, k), all consistently renormalizcd. Even 
though there is no explicit ghost content evident in the 
1-loop vertex Eq. , it does satisfy this identity at one- 
loop order p3 |. 

The dressed quark propagator appearing in Eq. is 
the solution to the gap equation, or the quark Dyson- 
Schwinger equation, which is 



,A 

S- 1 (p) = Z 2 S 1 (p) +C F Z 1 / g 2 D flv {p - q) 

J a 



7/i 



x S(q)r v (q,p), (8) 

where Sq (p) = ij ■ p + mb m , m^m is the bare current 
quark mass, and Z2(n 2 ,A 2 ) is the quark wave function 
renormalization constant. The general form for S(p)~ 1 
is 



Sip)- 1 =i~fpA(p 2 ,p 2 ) + B(p 2 ^ 2 ) 



(9) 



fi 2 is 



and the renormalization condition at scale p 2 
Sip) -1 *7 • P + rn(p,) where m(/z) is the renormalized 
current quark mass. 

Prior to the recent appearance of quenched lattice- 
QCD data pol l2l| . there had been little information 
available on the infrared structure of the gluon-quark ver- 
tex. The two 0{g 2 ) diagrams of Eq. @ can not be ex- 
pected to be adequate there. A common assumption [l6| 



has been to adopt an Abelian vertex Ansatz, such as 
the Ball-Chiu ^3 or Curtis-Pennington [l^ forms and 
attach the appropriate color matrix. In the case of an 
Abelian U(l) gauge theory, the counterpart to Eq. (JJJ is 
the Ward-Takahashi identity (WTI) 



k^,if fJ ,{p + k,p) = S(p + k)~ 



s( P y 



(10) 



At k = 0, the Abelian vertex f M has the same gen- 
eral form as given earlier in Eq. QJ. The Ward iden- 
tity t a {p) = -idS^ 1 {p)/dp a yields: &i = A(p 2 ), a 2 = 
2A'(p 2 ), and a 3 = 2B'(p 2 ), where /' = df{p 2 )/dp 2 . 
However, even if the Abelian Ansatz, igt c T a {p) 1 were 
to be adopted for the gluon vertex, it would not help in 
the present context because we need a representation in 
terms of an explicit set of Feynman diagrams for the re- 
sulting self-energy, in order to determine the symmetry- 
preserving BSE kernel. 

In Ref. [12j a study was made of a ladder summation 
Ansatz for the gluon vertex based on just the Abelian- 
like gluon exchange diagram of Eq. J2J; the symmetry- 
preserving BSE kernel was generated and used to explore 
meson and diquark masses. The vertex was generated by 
iterative and recursive techniques and, after convergence, 
is equivalent to solution of the integral equation 

+ k,p) = 2i7a - (C F - -yO J g 2 D^ v (p- q)7 M 
xS(q + k)T a (q + k,q)S(q) ll/ . (11) 

Here, at any order of iteration, the quark propagator is 
calculated by using the same vertex in the gap equa- 
tion, Eq. JHJ. Is this ladder sum a good approximation 
to the gluon-quark vertex, particularly in the infrared? 
The quenched lattice-QCD data indicates that the an- 
swer is no. The lattice data clearly gives a\{p 2 ) > 1 
for all available p 2 and the infrared limit appears to be 
ai(0) > 2.2. The ladder summation based on Eq. I|ll|l 
gives ot\{p 2 ) < 1, with infrared limit ai(0) ~ 0.94. The 
1-loop QCD analysis indicates that in the ultra-violet 
ai(p 2 ) approaches unity from above |2.'"{| . while the re- 
cent model vertex [22j, based on a non-perturbative ex- 
tension of the two 1-loop diagrams from Eq. (01, yields 
a\(p 2 ) > 1 for all p 2 , and agrees quite well with the lat- 
tice data. 

The reason for this problem can be seen from the 
color factors associated with the two 1-loop diagrams, 
Eq. © and Eq. Q), which are the leading terms in the 
ultra-violet region. The ladder sum in Eq. HTJ is built 
on the least significant of the two diagrams; the color 
factor of the omitted 3-gluon term is — N 2 times that 
of the retained term. The relative contribution to the 
Slavnov- Taylor identity, Eq. Q , from that term is of the 
same order at 1-loop. More generally, as discussed in 
Ref. [3, if G(fc 2 )(l - B(p, k)) > persists into the non- 
perturbative region, one can expect a\(p 2 ) > 1. One can 
also expect to obtain the wrong sign for a\{p 2 ) — 1 if a 
model kernel has the wrong sign. This is the case with 
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the Abelian-like ladder sum, Eq. Qllfl. Note that in an 
Abelian U(l) gauge theory, e.g., the photon-quark ver- 
tex, ai(p 2 ) = A(p 2 ) > 1. An Abelian Ansatz for this 
amplitude of the gluon-quark vertex might be quite rea- 
sonable, but it cannot be simulated by an explicit ladder 
sum — the color algebra prevents it. In analogy with the 
photon-quark vertex, where ai(p 2 ) > 1 is correlated with 
the spectral density being positive definite as the time- 
like region is approached, the gluon-quark vertex dress- 
ing has been referred to as an attractive effect in the 
infrared spacelike region 01 . (Of course, for the gluon 
vertex there should be no color octet bound states and 
no positive spectral density in the timelike region.) The 
3-gluon coupling is a strong source of the attraction at 
low spacelike p 2 ; it is N 2 times larger than the small 
repulsive effect of gluon exchange. 

The model for that we employ in this work, de- 
scribed in Section III, allows us to focus on zero gluon 
momentum. In this discussed and utilized in 

Ref. [22|, the two pQCD 1-loop diagrams for the ver- 
tex, Eq. 10 and Eq. 10} , are both closely related to the 
momentum derivative of the corresponding quark self- 
energy, apart from the differing color factors. The result- 
ing dependence upon the single quark momentum vari- 
able is similar for each diagram. Both are 1-loop integrals 
projected onto the same Dirac structures. We adopt the 
approach of Ref. to the vertex for our algebraic study; 
the approach is defined by taking the momentum depen- 
dence to be similar even in the infrared and with dressed 
propagators. Thus we combine the two terms and write 
Eq. J2| as 



rW(p+k,p) « Z lla -CC F / g 2 D^(p-q)^S (q+k) 

Jq 

X7^o(?)7„ (12) 



alent to the integral equation 

T a (p + k,p) = Zi-y a -CC F I g 2 D^(p~q)j fJi S{q + k) 

Jq 

xT a (q + k,q)S(q)j„, (13) 

is a natural suggestion. This was studied in Ref. [l4| . 
with S(q) being the self-consistent solution of the quark 
DSE, Eq. (23), containing the same dressed vertex. A fit 
to the lattice-QCD data for the vertex gave C = 0.51, 
a value that confirms that attraction by a mechanism 
outside the scope of iterated gluon exchange is present. 

An iterative representation is useful: = Yli=o 
where = Zij^, and i labels the number of internal 
gluon lines. The contribution with i + 1 internal gluon 
lines is obtained from the i th contribution by adding one 
gluon ladder. This is schematically depicted in Fig. ^ 




FIG. 1: (Color Online) The iterative relation for successive 
terms in the ladder- summed vertex. Here the large filled cir- 
cles denote the dressed quark-gluon vertex, the numbers in 
the parenthesis denote the numbers of gluon lines contribut- 
ing to the particular vertices and the small filled circles denote 
that the propagators are fully dressed. Note that an impor- 
tant non-Abelian term is approximately accounted for by the 
effective color factor C as described in the text. 



with C being an effective color factor to be determined by 
a fit to lattice-QCD data for the vertex. If the momentum 
dependence of the two combined terms from Eq. @ was 
identical, then we see that C = 1; this is equivalent to 
the Abelian limit. If one were to omit the 3-g luon term 
altogether, as in the iterative study in Ref. then C = 
(C F - ^f) Cp -1 , which for N c = 3, gives C = -1/8. One 
expects that the non-Abelian term is necessary for an 
effective model and thus that < C < 1. 

This vertex Ansatz allows us to avoid making a model 
for the dressed 3-gluon vertex for which there is little 
in the way of reliable information. It is implicitly hoped 
that the fit of C to lattice data will effectively compensate 
for deficiencies. Our aim is not the vertex itself but a 
study of the relative importance of classes of diagrams 
for the BSE kernel for meson masses. This vertex Ansatz 
allows an algebraic approach to the BSE meson masses 
that is quite illustrative of new qualitative information. 

From Eq. (|12f> . the non-perturbative summation equiv- 



A. A wider class of quark-gluon vertex dressing 

The enlarged class of dressing diagrams considered in 
this work is obtained iteratively as depicted in Fig. [21 
The contribution with i internal gluon lines is gener- 
ated from three contributions having a smaller number 
of gluon lines by adding one gluon ladder with dressed 
vertices. If the number of gluon lines in the three vertex 
contributions are denoted j, k and I, then summation is 
made over j, k and I such that j + k + I + 1 = i. Again, 
r M = Yli=o The iterative scheme is described by 

,A 

Ti(p + k,p) = -CC F ]T / g 2 D rJV {p-q) 

j.k,i J i 
i=j+k+i+i 

x Ti(p + k,q + k)S(q + k)Tl(q + k, q)S(q)T k v (q,p), 

(14) 
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for i > 1. 

If the iteration is carried to all orders, the equivalent 
integral equation is 

T^{p+k,p) = Zi7 M -CC F / g 2 D au {p-q)T a {p+k,q+k) 

Jq 

X S(q + k)T4q + k, q)S(q)T u (q,p). (15) 

If the iteration is stopped to produce all vertex functions 
with up to n internal two point gluon lines, our improved 
scheme takes into account l + n(n + l)(n + 2)/6 diagrams; 
the corresponding ladder-summed vertex at that order 
contains a subset of (n + 1) of these diagrams. 




FIG. 2: (Color Online) The iterative relation for the enlarged 
class of dressing diagrams considered in this work. Here the 
large filled circles denote the dressed quark-gluon vertex, the 
numbers in the parenthesis denote the numbers of gluon lines 
contributing to the particular vertices (with j + k + 1 + 1 = i) 
and the small filled circles denote that the propagators are 
fully dressed. The vertex contribution with i internal gluon 
lines is obtained from vertex contributions with less gluon 
lines. Note that an important effect of the non-Abelian 3- 
gluon coupling is approximately accounted for by the effective 
color factor C as described in the text. 



In Fig. |3 we use low order diagrams to illustrate the 
more general class of dressing terms included this way. 
Note that the included diagrams are restricted to planar 
diagrams. The contribution from crossed gluon lines in 
Fig. is not included. All diagrams of the ladder sum 
used in Ref. [l4|, such as Fig. E^, are included; the new 
element here is the self-consistent dressing of the internal 
vertices illustrated by Fig. |3b and Fig. Et. 

B. Symmetry-preserving Bethe-Salpeter kernel 

The renormalized homogeneous Bethe-Salpeter equa- 
tion (BSE) for the quark-antiquark channel, denoted by 
M, can be compactly expressed as 

[^M(k;P)] EF = / [K(k,q;P)]^[ X M(q;P)] GH , 

Jq 

(16) 

where rjf(fc;P) is the meson Bethe-Salpeter ampli- 
tude (BSA), k is the relative momentum of the quark- 
antiquark pair and P is their total momentum; E,...,H 




FIG. 3: (Color Online) Vertex skeleton diagrams at 0(g 5 ). 
Here the large filled circle denotes the quark-gluon vertex 
function dressed to the order of two effective gluon kernel 
lines. The small filled circles denote that the propagators are 
fully dressed. Previous work included the ladder structure 
typified by part (a). The enlarged class of dressing diagrams 
implemented in this work includes parts (b) and (c) as well. 
Non-planar diagrams such as part (d) are not accommodated 
by the present approach. We use an effective color factor to 
accommodate a major non-Abelian effect from the 3-gluon 
coupling as described in the text. 



represent color, flavor and spinor indices and the BS 
wavefunction is 

XM(k;P) = S(k+)r M (k;P)S{k-), (17) 

where k± = k ± and K is the amputated quark- 
antiquark scattering kernel. In general the kernel K is 
given in terms of the quark self-energy, S, by a functional 
relation dictated by chiral symmetry [lfj . This preserves 
the Ward-Takahashi identity for the color singlet axial 
vector vertex and ensures that chiral pseudoscalars will 
remain massless, independent of model details. 

In a flavor non-singlet channel, and with equal mass 
quarks, the axial-vector Ward-Takahashi identity is 

- ip^I ( P + p,p) = s-\p + p) 75 + Tsfi^GO 

-2m(n)r 5 (p + P,p), (18) 

where we have factored out the explicit flavor matrix. 
The color-singlet quantities and T 5 are the axial- 
vector vertex and the pseudoscalar vertex, respectively, 
and P is the total momentum. The amplitude T^(p + 
P,p) has a pseudoscalar meson pole. A consequence is 
that the meson BSE (|16f) for the (dominant) 75 amplitude 
at P 2 — is equivalent to the chiral limit quark DSE for 
B(p 2 ) and a non-zero value for the latter (DCSB) neces- 
sarily produces a massless pseudoscalar bound state |7J. 

The general relation between the BSE kernel, K, and 
the quark self-energy, S, can be expressed through the 
functional derivative |10| 

K(x', y'; x, y) = - J E(x', y'). (19) 

It is to be understood that this procedure is defined in 
the presence of a bilocal external source for qq and thus S 
and £ are not translationally invariant until the source 
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is set to zero after the differentiation. An appropriate 
formulation is the Cornwall- Jackiw-Tomboulis effective 
action [24|]. In this context, the above coordinate space 
formulation ensures the correct number of independent 
space-time variables will be manifest. Fourier transfor- 
mation of that 4-point function to momentum represen- 
tation produces K(jo, q; P) having the correct momentum 
flow appropriate to the BSE kernel for total momentum 
P. 

The constructive scheme of Ref. 0] is an example of 
this relation as applied order by order to a Feynman di- 
agram expansion for £(p). An internal quark propaga- 
tor S(q) is removed and the momentum flow is adjusted 
to account for injection of momentum P at that point. 
The number of such contributions coming from one self- 
energy diagram is the number of internal quark propa- 
gators. Hence the rainbow self-energy generates the lad- 
der BSE kernel. A 2-loop self-energy diagram (i.e., from 
1-loop vertex dressing) generates 3 terms for the BSE 
kernel. One can confirm that the axial-vector Ward- 
Takahashi identity is preserved. Similarly, the vector 
Ward-Takahashi identity is also preserved. 

To be more specific, with the discrete indices made 
explicit, we apply 



K 



GH 
EF 



5T, 



EF 



SS GH 



(20) 



to the self-energy given by the second term on the RHS 
of Eq. (JSJ • After a decomposition 



E(fe) = £s n (fc), 



(21) 



according to the number n of gluon kernels in the vertex 
defined by 

,A 

£"(fc)=C F / g 2 D p „{k-q) lp S{q)T^k), (22) 

Jq 

for n > 1, with 

,A 

S°(fc) = m bm + C F / g 2 D flv (k - q)l^S{q) lv , (23) 



The order n contribution to the BSE kernel is 

[K n (k,q;P)]%% = 

- C F g 2 D pv {k - q) [ 7 „] BG [r™(<Z-, k-)] HF 

-Cf^ g 2 D pi/ (k-l)[ lp S(l + )] EL 

x ^ *, -A m-,k.)] LF . (24) 

This format is the same as used in Refs. ^3 and [T^| . 
except that here the content of T™ is more extensive. 
With a bare vertex the first term of Eq. I|24|) produces the 
ladder kernel and the second term is zero. With a vertex 



up to 1-loop (n = 1), the first term of Eq. (|24|l produces 
the ladder term plus a 1-loop correction to one vertex; the 
second term produces two terms: a 1-loop correction to 
the other vertex and a non-planar term corresponding to 
crossed gluon lines. These three corrections to the ladder 
kernel have the same structure as the kernels shown in 
parts (b), (c), and (d) of Fig. El At higher order, n > 1, 
the BSE kernel produced in the present work departs 
from that considered in Ref. 

After substitution of Eq. into the BSE Eq. JTBJl, 
and with a change of variables, the meson BSE becomes 



,A 

T M (k;P) = -C F / g 2 D^(k-q)^ 

J a 



x \x M (q;P)Tu(q-,k-) + S(q+)A Mv {q,k;P)], (25) 

where we denote by Kmv the summation to all orders 
of the functional derivative of the vertex as indicated in 
Eq. lt24*|) . In particular, 



A M u(q,k;P) = J2A n Ml/ (Q,k;P), 



(26) 



with 



[^M,(q,k;P)} LF = 



K(q-,k-)} 



LF 



i SS G h(1±) 

x[xm(1;P)] gh . (27) 

The vertex iteration given in Eq. I|14|) produces the 
recurrence formula for the A^ 

Ab v (q,k;P) = -CC F ]T 

j,k,h 

J t A g 2 D p<J (q - t)Ti p (q + ,t + ) X M{t;P)T k u {t_,t- +k-q) 
xS(t- +k- q)T%(t- + k-q, kJ) 

+ j t A g 2 D p<T (k - t)r* p (q + ,t+ + q - k)S(t+ +q-k) 
xr*(t+ +q-k, t+) XM (t\ P)T*(t_, fc_) 

+ f t A g 2 D pa (q - t)Al Ip (q,t;P)S(t-)T k u (t-,t- + k - q) 
xS(t- + k- q)T*(t- +k-q, fc_) 

+ It" 9 2 D P Aq - t)Ti p (q + ,t + )S(t+)A k Mv {t, t + k-q-P) 
xS(t- +k- q)T*{t- +k-q, fc_) 

+ J t A 9 2 D pa {q - t)P p (q + ,t + )S(t + )T k (t + ,t + + k-q) 
xS(t+ + k- q)A h Ma (t + k-q, k; P)] 

(28) 

where A° Ml/ (q, k; P) = 0. 

The structure of the qq BS kernel produced by Eq. l|2*5^l 
and Eq. is schematically depicted in Figs. 0]and|S] 
With a general interaction kernel, g 2 D pa , it is exceed- 
ingly difficult to implement this formal recurrence rela- 
tion to obtain a BS kernel because of overlapping multiple 
integrals that compound rapidly with increasing order. 
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FIG. 4: (Color Online) Kernel decomposition. The filled tri- 
angles represent the meson BSAs, the filled circle represents 
the dressed quark-gluon vertex and the crossed circle repre- 
sents the A function. 




FIG. 5: (Color Online) A function decomposition. The filled 
triangles represent the meson BSAs, the filled circles repre- 
sent the dressed quark-gluon vertices and the crossed circles 
represent the A functions. The numbers in the parenthesis de- 
note the numbers of gluon lines contributing to the particular 
functions. 



III. ALGEBRAIC ANALYSIS 

A. The interaction model 

In the ultraviolet, the kernel of the quark DSE, Eq. JHJ), 
takes the form 

Z^D^T^p) -> 4ira(k 2 )^D f ™ c (kh„, (29) 

where k = p — q, and a(k 2 ) is the renormalized strong 
running coupling, which has absorbed the renormaliza- 
tion constants of the quark and gluon propagators and 
the vertex. The ladder-rainbow truncations that have 
been phenomenologically successful in recent years for 
light quark hadrons adopt the form of Eq. I|29|) for all k 2 
by replacing a(k 2 ) by a c s(k 2 ), which contains the cor- 
rect 1-loop QCD ultra-violet form and a parameterized 
infrared behavior fitted to one or more chiral observables 
such as (qq)^- In this sense, such an a c a(k 2 ) contains 
those infrared effects of the dressed vertex T„(q,p) that 
can be mapped into a single effective amplitude corre- 
sponding to 7„ for chiral quarks. Such a kernel does not 
have the explicit dependence upon quark mass that would 
occur if the vertex dressing were to be generated by an 
explicit Feynman diagram structure. In particular, one 
expects the vertex dressing to decrease with increasing 
quark mass; the effective ladder-rainbow kernel appro- 
priate to heavy quark hadrons should have less infrared 



strength from dressing than is the case for light quark 
hadrons. 

We use an explicit (but approximate) diagram- 
matic description of the dressed vertex T u (q,p), and 
to facilitate the analysis we make the replacement 
ATra oS (k 2 )/k 2 -> {2irfg 2 S 4 (k). This is the Munczek- 
Nemirovsky Ansatz 15] for the interaction kernel. The 
parameter Q 2 is a measure of the integrated kernel 
strength, and we expect this to be less than what would 
be necessary in ladder-rainbow format because of the in- 
frared structure now to be provided explicitly by the 
model vertex T )l (q 1 p). The equations of the previous 
Sections convert to model form by the replacement 

9 2 D^k) - (V - ^) (27r)W(*), (30) 

where we choose Landau gauge. It is the combination 
of Eq. (|30fl and the model vertex that is the DSE ker- 
nel; comparisons of Eq. (|30|) with information about the 
dressed gluon 2-point function are incomplete. The re- 
sulting DSEs for the quark propagator and gluon-quark 
vertex are ultra-violet finite; thus the renormalization 
constants are unity: Z\ = Zi = 1, and there is no distinc- 
tion between bare and renormalized quark current mass. 
We set rribm, = m(p,) = m. 

B. The algebraic vertex and quark propagator 

With this kernel, the vertex integral equation Eq. 115|) 
determines solutions for k — and we define T fl (p,p) := 
T^p). The resulting algebraic form for Eq. (|15|) is 

r„(p) = 7m - Cg 2 T a {p)S(p)Y ll {p)S{p)T a (p) . (31) 

In obtaining this form, we have used 3Cf/4 = 1, where 
the extra factor of 3/4 arises from the transverse projec- 
tor. The general form of the vertex is: 

r M 0) = ai(p 2 )7 A1 + a 2 {p 2 )l -pp^ -a 3 {p 2 )ip^ 

+ a^p 2 )^^ j-p (32) 

where a.i(p 2 ) are invariant amplitudes. From Eq. (|3 II) we 
find ola = 0, as was the case for the related models in 
Refs. 12] and [3. 

The vertex is a sum over contributions with exactly n 
internal effective gluon kernels according to 

oo 

r» = ]Tr», (33) 

71 = 

with the general contribution given by the recursive re- 
lation 

r» = -eg 2 m P )s(p)rl(p)s(p)rl( P ), 

j,k,l 
n=j + k+l+l 

(34) 
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FIG. 6: (Color Online) Quark propagator amplitude A(s) ver- 
sus Euclidean s — p 2 . We use the interaction mass scale Q = 
1 GeV and the current mass is m = 0.0183 Q = 18.3 MeV. 
C dependence calculated with converged summation of vertex 
dressing, for C — 0.15 (solid curve), C = 0.5 (dashed curve), 
C = —0.125 (dot-dashed curve) and C = — 0.25(dotted curve). 



where r?(p) = 7^. Substitution of the form S^p) -1 = 
17 • pA(p 2 ) + B(p 2 ) into Eq. (|3~4"|) gives T™(p 2 ) in terms of 
the functions A(p 2 ) and B(p 2 ). These latter functions 
must be solved simultaneously with the vertex at the 
given order. The algebraic form of the gap equation for 
the propagator is 



S-\p) 



ij ■p + m + g 2 ~/ fl S{p)T^(p), (35) 



where again the transverse projector and the color factor 
combine to yield 3Cf/4 = 1. After projection onto the 
two Dirac amplitudes we have 



A(p 2 ) = l-G 2 -tr 



B(p 2 )=m + g 2 -tr [7 p S(p)r»] 



(36) 
(37) 



Equations (|34[1 . i|36|) and (|37|l are solved simultaneously 
at a specified order n of vertex dressing. 

In the case where one is limited to a strict ladder sum- 
mation for the vertex with bare internal vertices, closed 
form expressions for the vertex amplitudes ctj in terms of 
A and B are obtainable 0, Q . With the enlarged class 
of dressing considered here, corresponding closed form 
expressions have not been obtained. However numerical 
evaluation is sufficient for the vertex and propagator am- 
plitudes; a numerical treatment of the BSE kernel must 



FIG. 7: (Color Online) Quark propagator amplitude A(s) ver- 
sus Euclidean s = p 2 . We use the interaction mass scale Q = 
1 GeV and the current mass is m = 0.0183 Q = 18.3 MeV. 
We show the influence of vertex dressing to order n as de- 
scribed in the text. For C = 0.15, n = yields the solid curve 
and the result is the ladder-rainbow truncation. The other 
curves are n = 1 (long dashed curve, 1-loop vertex), n = 2 
(dotted curve, 2-loop vertex), n = 12 (dot - short dashed 
curve), n = 13 (dot - long dashed curve) and n = 14 (dot 
dot dot - dashed curve) order of dressing of the quark gluon 
vertex. 



be made in any case. Numerical solution of the simul- 
taneous algebraic equations for the vertex and propaga- 
tor is carried out here using the algebraic and numerical 
tools of Mathematica (2{| with the assistance of the Feyn- 
Calc package used for computer-algebraic evaluation of 
the Dirac algebra pp| . 

The model parameter C for the vertex is determined by 
a fit to selected global features of quenched lattice-QCD 
data for the quark propagator |27j and the quark-gluon 
vertex 20]. This data is available for both quantities at 
current quark mass m — fh = 60 MeV. This is the same 
data as used to fit the same parameter C in Ref. [l4| : 
a different result will therefore reflect the wider class of 
vertex dressing herein. To facilitate comparison we also 
eliminate the role of the interaction strength mass scale 
parameter Q in this step by dealing with dimensionless 
quantities; Q will later be fixed by requiring that m p be 
reproduced. 

The lattice-QCD data for the quark propagator in- 
dicates that Z qu (0) = l/A qu (0) w 0.7 and M qu (0) = 
B(0)/A(0) w 0.42 GeV. Following Ref. [3 the lattice 
data for both the propagator and the vertex in the in- 
frared is characterized by the set of four dimensionless 
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FIG. 8: (Color Online) Quark mass function M(s) versus 
Euclidean s — p 2 . We use the interaction mass scale Q = 
1 GeV and the current mass is m = 0.0183 Q = 18.3 MeV. 
C dependence calculated with converged summation of vertex 
dressing, for C — 0.15 (solid curve), C = 0.5 (dashed curve), 
C = —0.125 (dot-dashed curve) and C = — 0.25(dotted curve). 



quantities evaluated at p 2 = 0: 

A(0,to 60 ) = 1.4 (38) 



FIG. 9: (Color Online) Quark mass function M(s) versus Eu- 
clidean s — p 2 . We use the interaction mass scale Q — 1 GeV 
and the current mass is m = 0.0183 Q = 18.3 MeV. We show 
the influence of vertex dressing to order n as described in the 
text. For C = 0.15, n = yields the solid curve and the result 
is the ladder-rainbow truncation. The other curves are n = 1 
(long dashed curve, 1-loop vertex), n — 2 (dotted curve, 2- 
loop vertex), n — 12 (dot - short dashed curve), n — 13 (dot 
- long dashed curve) and n — 14 (dot dot dot - dashed curve) 
order of dressing of quark gluon vertex. 



«i(0,m 60 ) = 2.1 (39) 



-Af(0,m 60 ) 2 a 2 (0,m 60 ) = 7.1 (40) 



-M(0,m 60 )a 3 (0,m 60 ) = 1.0, (41) 

where m 60 = fh/M qu (0). The best fit to these quan- 
tities gives C = 0.34 with an average relative error of 
f = 24 % and standard deviation ay = 70 %. The qual- 
ity of fit is about the same as in Ref. [l4|, an d changes 
AC « ±0.2 are not significant in this regard. For exam- 
ple, C = 0.15 leads to f = 39 % and oy = 72 %. We will 
use C = 0.15 because the resulting vertex at timelike p 2 is 
more convergent with respect to increasing order of dress- 
ing. The value of C being significantly greater than the 
strict ladder sum limit C = —1/8, we see that the attrac- 
tion provided by the 3-gluon coupling is important for 
the vertex. However the amount of attraction that must 
be provided in this phenomenological way in the present 
work is less than what was required in Ref. 14] to fit 



the same lattice quantities. In that work, C — 0.51 was 
found necessary. We attribute this difference to the fact 
that a wider class of self-consistent dressing diagrams is 
included in the present approach; attraction is provided 
by every vertex that is internal in the sense of Fig. [21 

In Figs. El and [3 we present the results for our calcula- 
tions of A(p 2 ) for different values of C and different orders 
of quark-gluon vertex dressing. We set Q = 1 GeV, so 
all dimensioned quantities are measured in units of Q. 
The current mass is m q = 0.0183 Q. One can see from 
Fig. © that C has a major impact on the behavior of 
A(s), especially in the timelike region. Fig.[7|shows that 
with n = 14 as the order of dressing of the quark-gluon 
vertex, we achieved convergence of the quark propagator 
function A(p 2 ) for p 2 > —Q 2 . The same is true for the 
function B(p 2 ). The relative measure of the convergence 
of the quark propagator functions with n is the conver- 
gence of the meson masses calculated using the solutions 
for the propagators. It will be shown later on that our 
calculations of and m p have converged to better than 
1% for n — 14. For heavier current quarks the conver- 
gence region for the solutions of A(p 2 ) and B(p 2 ) extends 
deeper into the time- like region of p 2 , which allows for 
convergent calculations of heavier meson masses. 
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In Figs. IS] and we present the results for our cal- 
culations of M(p 2 ) — B(p 2 )/A(p 2 ) for different values 
of C and different orders of quark-gluon vertex dressing. 
Again these calculations have Q = 1 GeV, so all dimen- 
sioned quantities are measured in units of Q. The vertex 
parameter C has a modest impact on the behavior of 
M(s). 

Figs. 1101 1111 and 1121 display the results for the vertex 
amplitudes aj(s) corresponding to different orders of ver- 
tex dressing. Successive orders after 1-loop (n = 1) serve 
to enhance the infrared strength for s < 1. The conver- 
gence with n is monotonic, in contrast to the convergence 
of the BSE kernel that is generated from this vertex, as 
discussed later. 

The quark condensate in the present model is given by 



3 



dss- 



B (s) 



Al(s)+Bi( S y 



(42) 



in terms of the chiral limit quark propagator amplitudes. 
There is no renormalization necessary because there is a 
spacelike sq for which Bq(s > sq) = 0. Because of the un- 
derrepresentation of the ultra-violet strength of the inter- 
action in this model, the condensate is characteristically 
too low. In particular we find 



-<W)c=o.i5 



(0.2146 Qf = (0.1266 GeV) 2 



(43) 



with Q = 0.59 GeV. The rainbow-ladder result (C = 0) 
is -(«9)lr = G 3 /{10n 2 ) = (0.1277 GeV) 3 . Thus one 
can see that the vertex dressing decreases the condensate 
slightly. In more detail, we have 



1.03, 



(44) 



:0.15 



which indicates that the ladder-rainbow truncation over- 
estimates the condensate by 3% compared to the more 
completely dressed vertex considered here. The previous 
study [l4| with a more restricted class of vertex dressing 
diagrams found that the ladder-rainbow truncation was 
18% too low. 



C. The algebraic Bethe-Salpeter kernel for mesons 

Substitution of the model interaction kernel Eq. 130(1 
into the meson BSE, Eq. I|25|l . produces the algebraic 
form 



T M (k;P) = -g 2 lll { XM {k;P)T ll {k_) 

+S(k+)A Mtl (k;P)}. (45) 

The previous general recurrence relation Eq. (|28fl for the 
general term of A Mv — S^Lo ^Mv now nas the algebraic 
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FIG. 10: (Color Online) Gluon-quark vertex amplitude a 1 (s) 
versus Euclidean s = p 2 , for C = 0.15. We use the in- 
teraction mass scale 5 = 1 GeV and the current mass is 
m = 0.0183 Q = 18.3 MeV. n = (solid curve) results from 
the bare vertex and is the ladder-rainbow truncation. The 
other curves are n = 1 (short dashed curve, 1-loop vertex 
dressing), n — 2 (dotted curve), n = 11 (dot - short dashed 
curve), n = 12 (dot - long dashed curve), n — 13 (dot dot dot 
- dashed curve) and n = 14 (long dashed curve). 



form 



Al Iv (k;P) = -Cg 2 



j,k,h 
n=j + k+h+l 

[P p (k + ) XM (k;P)T k u (k^)S(k.)T h p (k.) 
+Tj ) (k + )S(k + )T k u (k)xM(k;P)T h p (k.) 

+A J Mp (fc;P)s(fc_)rr£(fc_)s(fc_)r£(fc_) 

+ri(k + )S(q + )A k Mu (k;P)S(k.)T h p (k.) 

+rj ) (k + )S(k + )r k v (k + )S(k + )A h Ma (k ] p)} . (46) 

If we work at a given order, n, of vertex dressing, then 
the quark propagator, the dressed vertex, and the BSE 
kernel can be constructed recursively. By construction, 
chiral symmetry is preserved and the chiral pseudoscalar 
states are massless, independent of model parameters. 
Because of the algebraic structure, in which the BS am- 
plitude rjvf(fc; P) appears on both sides of Eq. (|45|l with 
the same qq relative momentum k, a physical solution 
where P 2 = —M^ is independent of k is defined only at 
k = 0. That is, the quark and antiquark have momenta 
r]P and (1 — rf)P. (Here we consider only equal mass 
quarks and thus have chosen r) = 1/2.) The Munczek- 
Nemirovsky model interaction does not allow momen- 
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FIG. 11: (Color Online) Gluon-quark vertex amplitude oi2(s) 
versus Euclidean s = p 2 , for C = 0.15. We use the in- 
teraction mass scale Q = 1 GeV and the current mass is 
m = 0.0183 = 18.3 MeV. n = (solid curve) results from 
the bare vertex and is the ladder-rainbow truncation. The 
other curves are n = 1 (short dashed curve, 1-loop vertex 
dressing), n = 2 (dotted curve), n — 11 (dot - short dashed 
curve), n = 12 (dot - long dashed curve), n = 13 (dot dot dot 
- dashed curve) and n = 14 (long dashed curve) . 



FIG. 12: (Color Online) Gluon-quark vertex amplitude 0:3 (s) 
versus Euclidean s = p 2 , for C = 0.15. We use the in- 
teraction mass scale 5 = 1 GeV and the current mass is 
m = 0.0183 Q = 18.3 MeV. n = (solid curve) results from 
the bare vertex and is the ladder-rainbow truncation. The 
other curves are n = 1 (short dashed curve, 1-loop vertex 
dressing), n — 2 (dotted curve), n = 11 (dot - short dashed 
curve), n = 12 (dot - long dashed curve), n — 13 (dot dot dot 
- dashed curve) and n = 14 (long dashed curve). 



turn transfer to quarks. This is a restriction present in 
all hadron studies made within this model. We define 
Tm{P) = Tm(& = 0;P), after which the form in which 
we solve the BSE is 

t m {p) = -Q 2 i,s{^) |r M (p)5(-|)r M (-|) 

+A Mm (0;P)}. (47) 

IV. MESON MASSES AND RESULTS 

The general form of a meson BS amplitude can be writ- 
ten as 

T M (fc; P) = ]T 1C\k- P) P M {k\k ■ P; P 2 ), (48) 

i 

where the K?(fe;P) are a complete set of independent 
covariants constructed from Dirac matrices and mo- 
menta that transform in a manner specified by the quan- 
tum numbers of the meson under consideration. The 
p M (k 2 ,k ■ P;P 2 ) are the corresponding invariant am- 
plitudes. (We do not show explicitly the color singlet 
unit matrix.) The model BSE under consideration here, 



Eq. (|47|l . has relative momentum k = 0, and the set of 
covariants is reduced considerably. We have 

N 

r M (P) = ]Tr(p)/i f (p 2 ) (49) 

i=l 

and it is convenient to develop a set of projection opera- 
tors Vj that allow us to isolate each amplitude according 
to 

f M = Tr D [PjTm] ■ (50) 

Then projection of the BSE, Eq. (|47|l . yields the eigen- 
value equation 

/(P 2 )=H(P 2 )/(P 2 ), (51) 

where / = (/L, / 2 f • ■ • ) is a vector of invariant ampli- 
tudes and the matrix H(P 2 ) is an N x N representation 
of the kernel. 

The mass, Mm, of the lowest bound state is obtained 
from the highest negative value of P 2 for which 

det[H(P 2 )-/] p2+M , f=0 = 0. (52) 

This method, namely the solution of the characteristic 
polynomial for Eq. 1511) , has also been followed in earlier 
work of this type [121 1 1 1| . 
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A. The Pion 

The general form of the ir Bethe-Salpeter amplitude 
requires four covariants and is 

T„(Jt; P) = 75 [i ft +-f-Pfl+i-k k-Pfi + a pv k p P v j% 

(53) 

in terms of amplitudes /^(fc 2 , k ■ P; P 2 ). We do not show 
flavor dependence since we treat we treat u-quarks and 
d-quarks the same in all other respects. In the present 
case only two covariants survive and we have 

r.(P)-75[^(P 2 )+7-^( J P 2 )]. (54) 
Convenient projection operators in this case are 

Vi = -\n, ^2 = ^7-^75. (55) 

B. The Rho 



TABLE I: Effect of quark-gluon vertex dressing to order n 
upon the masses of the ir and p mesons (in GeV) . The ladder- 
rainbow (LR) truncation corresponds to n = 0, one loop ver- 
tex dressing corresponds to n = 1, etc, while the full model 
result (converged to 3 significant figures) is labeled n = oo. 
Also displayed for m p is the mass error, Am p , and the rela- 
tive mass error, Am p /m p , of the LR truncation of the present 
model compared to a previous model 1 1 J I based on a limited 
class of vertex dressing diagrams. The mass scale parame- 
ter is Q = 0.59 GeV, the current mass of the u/d-quark is 
m = 0.0183 5 = 11 MeV, and C = 0.15. 



Vertex Dressing 




m p 


Am p 


Amp 
m p 


Amp [14] 


n = (LR) 


0.140 


0.850 


+0.074 


+0.095 


+0.295 


n = 1 (1-loop) 


0.135 


0.759 


-0.017 


-0.022 




n = 2 


0.135 


0.781 


+0.005 


+0.006 


+0.096 


n = 3 


0.135 


0.772 


-0.004 


-0.005 


N/A 


n = 4 


0.135 


0.778 


+0.002 


+0.003 


N/A 


n = oo (full model) 


0.135 


0.776 


0.0 


0.0 


0.0 



The general form of the p Bethe-Salpeter amplitude 
requires eight transverse covariants and corresponding 
amplitudes. Specific choices that have been found con- 
venient in earlier work are given in Refs. 0, |2^| . In the 
present case, the most general form is simply 

T P „(P) = (V - lvfl{P 2 )+V»vPvf P {P 2 )- 

(56) 

Again a unit color matrix is understood and we treat u- 
quarks and d-quarks as the same. Convenient projection 
operators that isolate the amplitudes are 

7>i = tU- P2 = -^<J^P,- (57) 



C. Vertex Dressing for Light Quarks 

There are a total of three parameters: C = 0.15, which 
has already been set by the quenched lattice data for the 
quark propagator and the gluon-quark vertex, while the 
experimental m n and m p are used to set the other two: 
the interaction mass scale, Q = 0.59 GeV, and the current 
mass for the u/d quark, m = 0.0183 Q = 11 MeV. The 
fully dressed vertex model is used in these determina- 
tions. In practice, we require convergence to 3 significant 
figures for the masses and this is achieved with a vertex 
dressed to order n = 14. Table [I] shows how the vertex 
dressing influences and m p . 

To confirm that our constructed BSE kernel preserves 
chiral symmetry, we verified that to any order of vertex 
dressing, and with m = 0, the chiral pion is massless 
to the numerical accuracy considered. The physical m v 
is not fixed perfectly by the symmetry but is almost so. 
The explicit symmetry breaking by the current mass is 
sufficient to determine m w for all orders of vertex dressing 



except for a few % error in the ladder-rainbow truncation 
(n — 0). Since the same behavior was observed in earlier 
work of this nature [l2L IbH ] , this result is quite model- 
independent. 

The response of m p to increasing order of vertex dress- 
ing shows that the ladder-rainbow truncation is missing 
74 MeV of attraction compared to the full model re- 
sult. The error decreases with each added order of vertex 
dressing. The relative error in the ladder-rainbow mass 
is 9.5% in the present vertex model, compared to 29.5% 
in the vertex model of Ref . 01 • Here each diagram for 
the dressed vertex has each of its internal vertices dressed 
in a self-consistent way. There are evidently both attrac- 
tive and repulsive contributions at the various orders that 
combine to add less net attraction to the ladder-rainbow 
truncation than what was found in Ref. FLU. 



D. Current Quark Mass Dependence 

One expects the influence of vertex dressing to de- 
crease with increasing quark mass because of the internal 
quark propagators in the vertex. Thus the LR truncation 
should become more accurate for mesons involving heav- 
ier quarks. It is useful to quantify this for the following 
reason. Phenomenological LR kernels are capable of 
incorporating many realistic features of QCD modeling 
and have been developed to provide efficient descriptions 
of light quark mesons, their elastic and transition form 
factors, and decay constants. A parameterized LR kernel 
that reproduces the experimental and m p has, by def- 
inition, absorbed the effective dressing of the vertex. The 
present work suggests that this is an amount of vertex 
attraction worth 9.5% of the vector meson mass. How- 
ever this phenomenological representation of the dressing 
does not have an explicit dependence upon quark mass 
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TABLE II: Error of the ladder-rainbow truncation for equal 
quark mass vector mesons in the u/d-, s-, c-, and b-quark 
regions, according to calculated mass and effective binding 
energies (in GeV). The ladder-rainbow (LR) truncation cor- 
responds to order n = in vertex dressing and the full model 
result corresponds to vertex dressing to all orders, n = oo, in 
this model. The mass scale parameter is Q — 0.59 GeV, and 
C = 0.15. 





ladder-rainbow 


full model 


LR % error 




n = 


n = oo 


this model 


441 


m„, d = 0.011 










m p 


0.850 


0.776 


9.5% 


30% 


B£ p 


0.346 


0.311 


11% 




m s — 0.165 










1714, 


1.08 


1.02 


6.0% 


21% 


B£,p 


0.350 


0.320 


9.0% 




m c = 1.35 










mj/^ 


3.11 


3.09 


0.3% 


3.5% 


BE j 14, 


0.260 


0.260 


0% 




m b = 4.64 










my 


9.46 


9.46 


0% 


0% 


B£r 


0.100 


0.100 


0% 





that would occur if the vertex dressing were to be gen- 
erated by an explicit Feynman diagram structure. One 
would expect such a phenomenological LR kernel to be 
progressively too attractive when applied to meson with 
progressively heavier quarks. 

The present model provides an opportunity to explore 
how much of the final meson mass result is attributable 
to vertex dressing and how this varies with quark mass. 
In Tabic ITU wc display results for the ground state vector 
mesons in the u/d-, s-, c-, and b-quark regions for both 
rainbow-ladder truncation and the full model. Again 
the quark current masses are determined so that the full 
model reproduces experiment. We see that the amount 
by which the LR masses are too large decreases steadily 
with increasing quark mass, as expected. The LR trunca- 
tion here is missing 6% of attraction for compared to 
21% in the restricted class of dressing diagram consider- 
ably previously ^| . The LR truncation is quite accurate 
for the cc and bb vector states, as expected. 

For the larger quark masses, the meson mass is dom- 
inated by the sum of the quark masses. We also ex- 
press the results in a form that has this large mass scale 
removed. For each state in Table |nj we display an ef- 
fective binding energy defined as B£ = 2M q (0) — my, 
where M q (0) is the quark mass function obtained from 
the DSE solution at p 2 = 0, and my is the meson mass. 
Thus Mg(0) is being used as a rough measure of the con- 
stituent quark mass. The use of a single p 2 point may 
well be an overestimate of constituent masses. Further- 
more, our fitted current quark masses are on the upper 
edge of what is usually quoted at a renormalization scale 
of \i = 2 GeV [23 . Such an overestimate would be am- 



TABLE III: The masses of the equal quark mass vector and 
pseudoscalar mesons in the u/d-, s-, c-, and b-quark regions, 
and the current quark masses required to reproduce the ex- 
perimental vector meson masses. All are in units of GeV. The 
values of m ria and m Vb are predictions. Experimentally |29||. 
m„ c = 2.9797 ± 0.00015 and m nb = 9.30 ± 0.03. The ficti- 
tious pseudoscalar 0J"s is included for comparison with other 
studies [3. 





= 0.011 


m s — 


0.165 


m c 


= 1.35 


mi, 


= 4.64 


m p 


= 0.776 


ra j, - 


= 1.02 




= 3.09 


mx(is) 


= 9.46 


BE P 


= 0.311 


BE 4, = 


0.320 


BEj/4, - 


= 0.260 


B£ T ~- 


= 0.100 


m-n 


= 0.135 




= 0.61 


m„ c 


= 2.97 


m Vb 


= 9.43 


BEn 


= 0.953 


B£ - = 


0.727 


B£ Vc = 


= 0.380 


B£ Vb - 


= 0.130 



plificd in the infrared region via a DSE solution for the 
quark propagator. Nevertheless, a relative comparison 
should be meaningful. Tabic [H] shows the dependence 
of B£ upon the current quark mass for the fully dressed 
model and the ladder-rainbow truncation. On this basis, 
the relative amount of overbinding of the LR truncation 
is consistent with its relative lack of attraction with re- 
spect to the mass results. 

In Table ITTT1 we display the full model results for both 
the vector and pseudoscalar qq states. The masses for 
rjc and rjb are predictions. In the c- and b-quark regions, 
these results are essentially the same as those of Ref . , 
because the differences in the employed model of vertex 
dressing become irrelevant when any dressing contribu- 
tion is suppressed by the large mass of propagators in- 
ternal to the vertex. The systematics of the mass depen- 
dence of hyperfine splitting that spans the c- and b-quark 
regions, here and in earlier work [l4|. strongly suggests 
that the experimental value j^jj, m nb = 9.30 ± 0.03, is 
too low. 



V. SUMMARY 

We have taken advantage of an algebraic model to en- 
large the class of diagrams for the quark-gluon dressed 
vertex that can be incorporated into the Bethe-Salpeter 
kernel, while allowing a practical application to the calcu- 
lation of meson masses. A given expansion of the vertex 
in diagrammatic form, produces a diagrammatic expan- 
sion of the quark self-energy, which in turn specifies a 
diagrammatic expansion of the BSE kernel if chiral sym- 
metry is to be respected. This procedure relieves the 
phenomenology of the task of reproducing Goldstone's 
theorem whenever parameters are changed - it is always 
obeyed in this approach and thus phenomenology can ad- 
dress itself to a more constrained task. The constraints 
are considerable: a realistic ladder-rainbow kernel fitted 
to (qq) |3j produces m p , and itik * to better than 5%. 
Such a phenomenological LR kernel for light mesons has 
absorbed vertex dressing but without the explicit m q de- 
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pendence associated with an explicit diagrammatic rep- 
resentation of the dressed gluon-quark vertex. In order 
to gain more information it is necessary to work with a 
model that can implement a summation of vertex dia- 
grams, turn that into a summation of diagrams for the 
chiral symmetry preserving BSE kernel, and allow a prac- 
tical solution of the meson BSE. 

To this end we use the Munczek-Nemirovsky Ansatz 
[l5j for the interaction kernel. We use an improved model 
for the quark-gluon dressed vertex wherein each diagram 
for the dressed vertex has each of its internal vertices 
dressed in a self-consistent way. This moves consider- 
ably beyond the ladder BSE structure [l4| for the ver- 
tex, in which vertices internal to the dressed vertex of 
interest are bare. In common with Ref. ^3], we also 
use an effective method, with one parameter (C = 0.15 
for this model), to accommodate the important non- 
Abelian effect of the 3-gluon coupling for the vertex. 
Quenched lattice-QCD data for the quark propagator 
and the quark-gluon vertex at zero gluon momentum 
fixed the parameter C, while m v and m p fixed the other 
two parameters via the fully dressed vertex results. 

The resulting model provides a laboratory within 
which the relevance of ladder-rainbow truncation (bare 
vertex) can be explored over a range of quark masses 
from u/d-quarks to b-quarks. The influence of the en- 
larged class of vertex dressing diagrams included in this 
work is seen to indicate that LR truncation is missing 
9.5% of attraction for m p , whereas the previous informa- 
tion from a smaller class of vertex dressing diagrams had 
LR missing 30% of attraction. The extra dressing dia- 
grams included here tend to provide some cancellation, 
making the LR truncation somewhat more accurate. As 
heavier qq mesons are considered, the amount of miss- 



ing attraction in the LR truncation decreases steadily, as 
does the influence of vertex dressing — it is less than 1% 
for the J /if) and T. 

The influence of the non-Abelian 3-gluon coupling is 
very significant. No attempt has been made to consider 
4-gluon coupling nor to consider non-planar gluon line di- 
agrams (e.g. crossed-box diagrams) for the vertex dress- 
ing. On the other hand, a limited class of non-planar 
gluon line diagrams for the meson BSE kernel, as gener- 
ated from the planar diagrams of the dressed vertex, are 
included. While the complex task of including both pla- 
nar and non-planar two-point gluon line diagrams for the 
vertex is currently underway, it is not yet known whether 
explicit 3- and 4-gluon couplings can be accommodated, 
even through the device of an effective color factor. 
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